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ABSTRACT 

We identify a sample of 74 high-redshift quasars (z > 3) with weak emission lines from the Fifth Data Release 
of the Sloan Digital Sky Survey and present infrared, optical, and radio observations of a subsample of four 
objects at z > 4. These weak emission-line quasars (WLQs) constitute a prominent tail of the Lya+N V equiv- 
alent width distribution, and we compare them to quasars with more typical emission-line properties and to 
low-redshift active galactic nuclei with weak/absent emission lines, namely BL Lac objects. We find that 
WLQs exhibit hot (T ^ 1000 K) thermal dust emission and have rest-frame 0.1-5 fim spectral energy distri- 
butions that are quite similar to those of normal quasars. The variability, polarization, and radio properties of 
WLQs are also different from those of BL Lacs, making continuum boosting by a relativistic jet an unlikely 
physical interpretation. The most probable scenario for WLQs involves broad-line region properties that are 
physically distinct from those of normal quasars. 
Subject headings: quasars : general - quasars : emission lines 



1. INTRODUCTION 

The optical and ultraviolet (UV) spectra of type 1 quasars 
are characterized by a blue power-law continuum and 
strong, broad permitted emission lines with FWHM = 1000- 
20,000 km s"'. These lines are thought to originate in pho- 
toionized gas located ^ 10'^"'** cm from the central black 
hole, with differential gas velocities produ cing Doppler- 
broadened line profiles (e.g.. iPetersoiil 119971: l&o Hk 1999|). 
This gas may be confined to individual clouds outside of the 
accretion disk or it may be related to th e accretion disk itself, 
perhaps as part of a disk wind (e.g., lEmmering et al.l [19921; 
iMurrav et al.|[T995h . 

Observationally, the strongest of these emission lines, in 
terms of both flux and equivalent width (EW), is Lya A1216. 
Results from UV spectra at z < 2 and optical spectra at z > 2 
(e.g., Schneider et al. 1991", Francis et al. 1993; Osmeretal. 
1994; Warren et al. 1994; Zheng et al. 1997; Brotherton et al. 
20011: IDietrich et all 120021) show that the rest-frame EW of 



LyofH-N V A1240 is typically 50-110 A. These results, how- 
ever, are based on small numbers of sources, with tens 
to hundreds of quasars per study. The Sloan Digital Sky 
Survey (SDSS, York et al. 2000) has substantially increased 
the available spectroscopic sample of high-redshift quasars; 
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more than 5000 quasars have been discovered at z > 3 
(Schneider et al. 2007), extending all the way to z = 6.42 
dFan et al.l i2003). The SDSS quasar selection algorithm at 
these redshifts ( Richards et al. 2002a) is based largely on the 
red colors produced by the Lyman break (Aiest = 912 A) and 
the onset of the Lya forest (Arest = 1216 A). It is sensitive 
to quasars with bright UV continu a, without a stro ng de- 
pendence on emission-line strength dFan et al.l 120011) . and a 
small fraction of high-redshift quasars has been found with 
very weak or absent Lya emission (e.g iFan et all [T999t 
lAnderson et ani^OOTtlCoUinge et al.ll200llFan et alJI200^ 

The first weak-emission line quasar (WLQ) fou nd at high 
redsh ift was SDSS J153259.96-003944.1 (z = 4.62. IFan et afl 
119991 hereafter SDSSJ1532). Its flat continuum and lack of 
emission lines suggested that it could be the highest-redshift 
BL Lac object ever found, but its paucity of radio flux, X- 
ray flux, and optical polarization were inconsistent with the 
BL Lac hypothesis. The spectra of classical BL Lacs lack 
emission lines because their continuum emission is relativis- 
tically boosted by a jet along the line of sight, and they are 
usually radio-loud. X- ray-loud, and highly polarized (e.g., 
lUrrv & Padovani|[l995l) . The nature of SDSSJ1532 was left 
as an open question. 

Subsequently, two additional WLQs, 

SDSS J130216.13h-003032.1(z = 4.47, hereafter SDSSJ1302) 
and SDSS J144231.72-H01 1055.2 (z = 4.51, hereafte r 
SDSSJ1442), were identified by lAnderson etlTI (1200 ll) . 
and the strong radio and X-ray em ission of the latter 
dSchneider et all 120051: iShemmer etall |2006) indicated that 
its continuum could be moderately beamed. In a sample of 
optically selected BL Lac candidates, Collinge et al. (20051) 
found seven sources at z > 2.7, includin g three radio sources 
and one X-ray source dSchneider et al.l 12003) . The highest- 
redshift WLQ, SDSS J133550.8 1-H353315.8 at z = 5.93, was 
discovered by Fan et alj d2006l) . who discussed whether it 
could be a strongly lensed galaxy, a BL Lac object, or a 
quasar with a very weak emission-line region. 

Further insight into the nature of WLQs has been 
gained from pointed X-ray observations with Chandra. 
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IVignali et all (1200 Ih p re sented observations of SDSSJ1532, 
while Ishemmer et alj (12006b presented deeper data for 
SDSSJ1532, as well as observations of SDSSJ1302, 
SDSSJ1442, and S PSS J140850.91+0205 22.7 (z = 4.01, here- 
after SDSSJ1408); Shem mer etaP (20Q6) discussed whether 
the weak emission lines could be due to continuum boo sting 
or a deficit of line-emitting gas. IShemmer et alJ (l2009l) pre- 
sented deeper data for SDSSJ1302 as well as observations of 
eight additional WLQs at z > 2.7; they also compared the op- 
tical, X-ray, and radio properties of WLQs to those of BL 
Lacs, and discussed whether WLQs could be extreme quasars 
with high accretion rates. They found that WLQs are weaker 
than BL Lacs at radio and X-ray wavelengths relative to the 
optical, but that there is no evidence of steep hard X-ray spec- 
tra characteristic of high accretion rates. 

This paper extends the study of WLQs to the full high- 
redshift (z > 3) sample that has been discovered by the SDSS, 
and uses a rich, multiwavelength data set for a subsample of 
four WLQs at z > 4 to test the available hypotheses regard- 
ing their nature. In Section|2] we measure EW(LyaH-N v) for 
each z > 3 quasar in the Fifth Data Release Quasar Catalog 
(ISchneider et alJl2007 ). select quasars with EW values > 3a 
below the mean of the distribution, and compare their prop- 
erties to those of thegeneral high-redshift quasar population. 
Sections |3] H] and [Sfthen present results from an extensive 
observational campaign focusing on the WLQs SDSSJ1302, 
SDSSJ1408, SDSSJ1442, and SDSSJ1532. These data in- 
clude Spitzer 3.6-24 /im mid-infrared (mid-IR) photometry 
and near-IR photometry in 7, H, and bands that are com- 
pared to the spectral energy distributions (SEDs) of typical 
quasars (Section |3]i; multiple epochs of optical spectroscopy 
and /-band photometry (in addition to the available SDSS 
data) that constrain the line and continuum variability of 
WLQs (Section |4]i; and optical imaging polarimetry along 
with Very Large Array (VLA) radio continuum observations 
at L and C bands that test for the presence of synchrotron 
emission from a relativistic jet (Section|5]l. We discuss the im- 
plications of our results and the conn ection to lower-redshift 
quasars wit h weak emission lines (e.g. .' McDowell et al.ll995l : 
iLeighly eta l. 2007a) in Section |6] and we summarize and 
conclude in Section|7] 

We assume a cosmology with f^A = 0.7, Q,m = 0.3, and Hq = 
70 km s~' Mpc"'. All quoted emission-line EWs are in the 
rest frame of the quasar. 

2. THE WEAK EMISSION-LINE QUASAR SAMPLE 

The SDSS Fifth Data Release Quasar Catalog in- 
cludes 77,429 object s and covers an area of 5740 deg^ 
(ISchneider et al.ll2007h . The 5374 quasars at z > 3'" are se- 
lected down to a magnitude limit of / = 20.2 dRichards et alj 
l2006ah . The SDSS spectra have wavelength coverage A = 
3800-9200 A and resolution R - 2000. In Figure[T] we show 
six spectra exhibiting a range of redshifts, emission-line EWs, 
continuum spectral slopes, and signal-to-noise ratios (S/Ns), 
including two WLQs. 

For each spectrum, we fit a power law" of the form 
f\=CxX^ to regions of the continuum that are generally 

Lya is detectable in SDSS spectra down to j ~ 2.2, but we focus our 
analysis on z > 3 w here the flux limit of the quasar survey is a magnitude 
fainter (e.g., [Richar ds et al. 2006a, see their Figure 5). The highest-redshift 
quasar in the DR5 catalog is at z = 5.41. 

' ' The power-law spectral slope of quasars is often expressed in terms of 
a, where oc u". We use this definition of a throughout the paper. Here, 
P = -a-2. 
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Fig. 1. — Spectra of six z > 3 quasars exhibiting a range of redshifts, 
LyQ+N V emission-hne EWs, continuum spectral slopes, and S/N. The wave- 
length corresponding to Arest = 1216 A is marked by a dotted line in each 
spectrum. The objects in the second and fifth panels are WLQs. 



uncontaminated by emission lines (1285-1295, 1315-1325, 
1340-1375, 1425-1470, 1680-1710, 1975-2050, and 2150- 
2250 A). We then integrate the flux above the power-law con- 
tinuum between Arest = 1160 A and Arest = 1290 A to deter- 
mine emission-line flux and EW. This region includes the 
Lya A1216, N V A 1240, and Si II A 1263 emission lines, so 
we measure the strength of a blended emission feature that is 
dominated by the Lya and N V components. The median S/N 
of the spectra integrated over this wavelength region is 140 
for line+continuum and 95 for the continuum. 

The EW(LyQ!-i-N v) measurements for all 5374 quasars are 
given in Table [T] The typical statistical uncertainty of each 
measurement is ^ 1-2 A. There is also a systematic uncer- 
tainty in the placement of the continuum level, which depends 
on the wavelength regions chosen for the continuum fit. Vi- 
sual inspections show that the fits accurately describe the data, 
and while the EWs of individual objects do change slightly 
(~ 10%) if different wavelength regions are used to fit the 
continuum, the aggregate results are quite similar An impor- 
tant caveat is that our method underestimates the intrinsic line 
strength whe n strong intervenin g absorption is present. Some 
studies (e.g.. [Dietrich et al"]|2002i) interpolate across Lya for- 
est absorption features and fit the intrinsic emission-line pro- 
file, but this requires an accurate determination of the ab- 
sorption systems along the line of sight to each quasar and 
is difficult without high-S/N data at higher spectral resolu- 
tion. The effects of Lya forest absorption are discussed fur- 
ther below and in Section 12.1.31 In addition, there is a break 
in the quasar continuum near Arest = 1 100 A (e.g., Zheng et aTl 
[T997; .Shang et al. 2005) . so our power-law continuum fit at 
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Lya + NV rest-frame EW [A] 

Fig. 2. — Distribution of rest-frame Lya+N V EWs for tlie sample of 3058 z > 3 non-BAL quasars described in Section[2] The dotted red line sliows the best-fit 
log-normal distribution, with a mean EW of 63.6 A and a icr range of 39.6-101.9 A. The fit is a reasonable description of the data above EW = 20 A, but there is 
a prominent tail toward low EW values. The dashed line and the hashed section of the histogram indicate our WLQ definition. 



Arest > 1285 A may slightly overestimate the continuum flux 
at Lya and therefore underestimate the strength of the emis- 
sion feature. 

We restrict our analysis to the 3979 sources that were iden- 
tifi ed as primary quasar t argets by the final algorithm given 
by [Richards et alJ (l2002a|). These objec t s ha ve the uniform 
selection flag set in the I Schneider et alJ (|2007|) catalog and 
constitute a statistical color-selected sample. We addition- 
ally exclude 915 broad-absorption line (BAL) quasars from 
the catalogs of Trump et al. (2006) and Gibson et al. ( 2009), 
as well as six objects flagged by visual inspection as having 
either BAL-like absorption or missing data in the Lya region. 
The final sample includes a total of 3058 quasars. 

The distribution of Lya-i-N V rest-frame EWs for this sam- 
ple is shown in Figure |2] The median and mean of the mea- 
sured values are 62.3 A and 67.6 A, respectively. The distri- 
bution is well described by a log-normal function with best- 
fit parameters (logEW) = 1.803 and cr(logEW) = 0.205, cor- 
responding to a mean of 63.6 A and a la range of 39.6- 
101.9 A. However, there is a prominent tail toward small 
values. If the log(EW) values were indeed distributed like 
a Gaussian, where the 3cr range encompasses 99.73% of the 
curve area, we would expect to find ~ 4 objects beyond this 



range (EW= 15.4-261.9 A) at each end. Instead we find none 
at the high end and 56 at the low end, and we define our sam- 
ple of WLQs to be these 56 objects with EW < 15.4 A. This 
WLQ definition is more inclusiv e than the EW < 10 A defini- 
tion adopted by Shemmer et al.l (12009). but the precise cutoff 
between WLQs and "normal" quasars is somewhat arbitrary; 
our goal is simply to select the sources with extreme EWs 
that exist beyond what would be expected from a simple log- 
normal distribution. 

To illustrate the spectroscopic properties of WLQs relative 
to the rest of the quasar sample, we construct composite spec- 
tra in bins of EW(Lya-i-N v) and present them in Figure [3] 
These composite spectra span the range Ar est = 800-1800 A 
and are constructed following the method of lFan et 311(120041) : 
each spectrum is normalized based on its flux at Arest = 1450 A 
and the average of all spectra is computed. We use five EW 
bins: (1) EW > la above the mean, 428 objects: (2) EW be- 
tween the mean and la above it, 1038 objects: (3) EW be- 
tween the mean and Icr below it, 1000 objects: (4) EW be- 
tween la and 3a below the mean, 536 objects: and (5) EW 
> 3a below the mean, 56 objects. In the first three compos- 
ite spectra, the red side of the Lya-i-N V blend (dominated 
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Fig. 3. — Composite spectra of z > 3 quasars in different EW(Lya+N V) 
bins. The wavelengths of Lya and C IV are marked by dotted lines. 
Top panel: Composite for 428 quasars with EW > Ict above the mean 
(EW > 101.9 A). Second panel: composite for 1038 quasars with EW be- 
tween the mean and Icr above it (63.6 A < EW < 101.9 A). Third panel: 
composite for 1000 quasars with EW between the mean and Icr below it 
(39.6 A < EW < 63.6 A). Fourth panel: composite for 536 quasars with EWs 
between la and 3cr below the mean (15.4 A < EW < 39.6 A). Fifth panel: 
composite for 56 WLQs with EW > 3(t below the mean (EW < 15.4 A). 
Bottom panel: composite for the subset of 32 WLQs that do not exhibit sig- 
nificant C IV emission or PDLA absorption (see Table|2). There is a trend for 
the Lyo/N V ratio to decrease toward lower EWs. 



by N v) is nearly identical, but the blue side (dominated by 
Lya) steadily decreases in strength; while in the fourth spec- 
trum, the peak flux shifts redward of Arest = 1216 A as the 
N V component begins to dominate. This trend continues in 
the WLQ composite, where the 1216 A flux actually drops 
below the (unabsorbed) continuum level. This extreme be- 
havior is driven by a handful of WLQs with strong absorp- 
tion features, and in f act 15/56 WLQs are included in the 
iProchaska et alj (l2008h catalog of proximate damped Lya ab- 
sorbers (PDLAs, Av < 3000 km s"', A^^ > 2 x 10^° cm'^). 
In the bottom panel, we exclude these WLQs whose weak 
Lya-i-N V lines can be explained by absorption rather than in- 
trinsic weakness, as well as 10 additional WLQs that exhibit 
EW(C IV) > 10 A (these sources ai-e flagged in TableEll. The 
resulting WLQ composite spectrum has flux roughly equal to 
the continuum level near Lya A1216, a weak emission feature 
near N V A1240, and a neghgible bump near C IV A1549. 

There is also a trend for the continuum luminosity to in- 
crease in the lower EW bins. The median continuum lumi- 
nosity for the whole sample is 41% larger than the median 
for sources in the highest EW bin and 39% smaller than the 
median for WLQs. This trend is consistent with the Baldwin 
effect (the observed anti-correlation between C iv EW and 
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Fig. 4.— Distiibution of C IV EWs for 2035/3058 quasars from Figure[2] 
measured from emission-line fits performed by Shen et al. ( 2008). The dotted 
red line shows the best-fit log-normal distribution, with a mean EW of 4 1 .9 A 
and a Ict range of 26.1-67.4 A. Objects with the weakest C IV lines are not 
included because accurate EWs could not be determined. 



quasar luminosity: iBaldwinll 1 9771) . 

We select the WLQ sample based on extreme 
Lya-i-N V EWs rather than C IV because the former line is 
stronger, and therefore readily detected even at low S/N, 
and observable at the highest redshifts (z > 4.9). However, 
we do compile IV EWs measured from emission-line fits 
performed by iShen et alj (2008). These values are presented 
in Table [T] and their distribution is shown in Figure |4] 
For the 2035/30 58 quasars in our final sample where the 
IShen et all (12008 ) fits yield secure C IV EWs, the best-fit 
log-normal distribution has parameters (logEW) = 1.622 and 
cr(logEW) = 0.206 (mean41.9 A, larange 26.1-67.4A), and 
the median EW ratio (Lya-i-N V)/C IV is 1.53. This ratio is 
smaller than what has b een reported i n prev i ous studies (e.g., 
'Schneider et al.' '1991'; 'Osmer et al.' '1994^; ' Brotherton elaLl 
2001; Vanden Berk et al. 2001; Dietrich et alj |2002|), which 
tend to find larger EW(Lya-i-N v) and smaller or comparable 
EW(C IV) values. As mentioned above, our method under- 
estimates the intrinsic strength of Lya-i-N V emi ssion when 
strong Lya forest absorption is present, and the IShen et all 
(2008) fits overestimate EW(C IV) when strong Fe II emission 
mimics a broad Gaussian component of the C IV profile; but 
our results for both lines are nonetheless within the range of 
litera ture values (e.g., Francis et al. 1991). Also, early studies 
(e.g.. I Schneider et"ar.i,1991; Osmer et al.1ll994h were biased 
toward more luminous quasars, so their lower EW(C IV) 
values could be explained by the Baldwin effect. There is 
no information about objects with the weakest C IV lines 
in Table [T] or Figure |4] b ecause accurate EWs could not be 
measured from the Sh en et al.l (12008) fits; in fact, 34/56 
WLQs and all objects with EW(C IV) < 6 A ai-e not included. 

2. 1 . Sample Properties 

To investigate the nature of the low-EW tail of the distribu- 
tion, in this section we compare several properties of WLQs 
to those of the general z > 3 quasar population. We perform 
this analysis on all 56 WLQs in the uniform sample and on the 
subset of 3 1 WLQs with no C IV or PDLA flags (see above 
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and Table |2|l; the results are insensitive to the choice of sam- 
ple, and we quote values for the full uniform sample. 

2.1.1. Emission-line and Continuum Properties 

In addition to EW, our continuum fits to the SDSS spectra 
yield values of LyaH-N V emission-line flux, continuum flux, 
and continuum spectral slope for each source. We find that 
the emission-line luminosities of WLQs are all at the low end 
of the distribution for normal'^ z > 3 quasars (28/56 WLQs 
are in the bottom 1%, 44/56 are in the bottom 5%, and all 
are in the bottom 26%); by the Kolmogorov-Smirnov (K-S) 
test, the probability that the two samples are drawn from the 
same parent distribution is /:> = 4 x 10"-^'. The median line 
luminosity for WLQs is a factor of 3.8 lower than the me- 
dian for normal quasars, indicating that Lya is weak in an 
absolute sense, rather than just relative to the continuum. As 
mentioned above, the continuum luminosities for WLQs are 
somewhat larger on average than those for normal quasars, but 
their continuum slopes are quite similar. We find median val- 
ues of a = -0.54 and a = -0.52 for normal quasars and WLQs, 
respectively, and the K-S test gives a probability p = 0.20 that 
the distributions are the same. 

2.1.2. Radio Flux 

We also investigate the radio properties of WLQs using 
data from the Faint Images of the Radio Sky at Tweny cm 
( FIR,ST) surv ey (Becker et al. 1995). A search of the catalog 
(1 White et al.lll997l) within 2" of each quasar yields 147 de- 
tections and 2809 upper limits'^; 102 quasars, including four 
WLQs, are not covered by the FIRST footprint. We find that 
the radio-detection fraction is 4.9% (142/2904) for normal 
quasars and 10% (5/52) for WLQs. This difference is not sig- 
nificant; given the radio-detection fraction of normal quasars, 
the binomial distribution gives a 1 1 % chance of having five 
or more detections in a sample of 52. If we relax the uniform 
selection criterion described above and include all non-BAL 
quasars, the radio fraction increases to 6.1% (233/3811) for 
the general population and to 21% (15/70) for WLQs. This 
increase, which is especially prominent for WLQs, is due to 
the fact that FIRST sources were targeted for spectroscopy by 
the SDSS even if their optical colors did not meet the final 
quasar-selection criteria. Assuming a 6.1% detection proba- 
bility, the chance of having 15 or more detections in a sample 
of 70 is only 0.002%. 

A more relevant quantity than the radio detection fraction 
is the radio-loudness parameter R, which is defined as the ra- 
tio of radio to optical flux density. We adopt the definition'"^ 
used bv LTiang et al.1 (IM)! . R = /^(6 cm)//^(2500 A). We 
calculate the rest-frame 6 cm flux density from the observed 
20 cm flux density assuming a power-law slope ar = -0.8, and 
we calculate the rest-frame 2500 A flux density from the ob- 
served z band (A = 893 1 A) flux density assuming a power-law 

We use the term "normal" hereafter to refer to non-WLQs (i.e., quasars 
withEW> 15.4 A). 

We calculate the detection threshold corresponding to the rms at each 
position. This is a 5fT upper limit, including CLEAN bias: Fy = 5(T+0.25 mJy 
(White et al.|[T997l) . 

The original definition of R by Keller marm et alj (19891 measures the 
optical flux density at the B band, A = 4400 A. We choose a shorter wave- 
length to more closely match the rest-frame wavelengths observed by the 
SDSS. For a power-law slope a = -0.5, one would expect 33% more flux 
density at 4400 A than at 2500 A so the definition we use generally produces 
larger R values. 
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redshift 

Fig. 5. — Redshift distributions of normal z > 3 quasars (solid black line, 
bin size Az = 0.05) and WLQs (dashed red line, bin size Az = 0.20). The two 
distributions diverge at z ~ 4.2, with WLQs appearing to be more common 
at higher redshift. The i nset panel shows the selection probability for each 
WLQ (see Section l2.1.3t . indicating that weak-lined sources are more likely 
to be selected at higher redshift. The increased opacity of the IGM near 
■^rest = 1216 A at higher redshift is also likely to bias EW measurements 
toward lower values. 



slope ao = -0.5. We find that 3.2% (92/2904) of normal z > 3 
quasars and 6% (3/52) of WLQs have R > 100, while > 4.6% 
(133/2904) and > 10% (5/52) have R > 10. The latter frac- 
tions are lower limits because the FIRST detection threshold, 
typically ~ 1 mJy, is only deep enough to constrain R <10 for 
the optically brightest sources (~ 15% of the sample). Again, 
including all non-BAL quasars changes the R > 10 fraction 
to > 5.9% (224/3811) and > 21% (15/70), and the R > 100 
fraction to 4.1% (156/381 1) and 7% (5/70). 

To compare the distributions of R values statistically, given 
that most values are upper limits (i.e., the distributions are 
heavily censored), we u se the survival analysi s software pack- 
age ASURV Rev 1.3 (iLavaflev et al.1 Il992h. which imple- 
ments the methods presented in iFeigelson & NelsonI ([T985h . 
We perform the Gehan, logrank, and Peto-Prentice two- 
sample tests. In the uniformly selected sample, these tests 
give probabilities p = 0.097, p = 0.145, and p = 0.147 that the 
WLQs and normal z > 3 quasars have the same distribution of 
radio loudness. If the uniform-selection criterion is dropped 
(i.e., all non-BAL quasars are included), these probabilities 
drop to p=lx 10-^ p = 2x 10-^ and = 8 x 10"^ Thus, 
there is no evidence of a difference in radio properties be- 
tween WLQs and normal quasars among the sources that were 
selected for spectroscopy solely on the basis of optical colors. 
However, the use of FIRST data to identify SDSS spectro- 
scopic targets selects a larger fraction of WLQs; of the non- 
BAL quasars at z > 3 with FIRST detections, 6.0% (15/248) 
have EW(Lya+N v) < 15.4 A. This may be because radio se- 
lection is not biased against quasars with weak Lya emiss ion, 
while color selection at z < 4 is biased (see Section l2.1.3b . 

2.1.3. Redshift Distribution 

The redshift distributions of normal z > 3 quasars and 
WLQs are shown in Figure |5] The two distributions di- 
verge at z ~ 4.2, with WLQs constituting 1.3% (36/2737) of 
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the quasar population below this redshift and 6.2% (20/321) 
above it. The redshifts measured for individual WLQs can 
be uncertain by a s m uch as & = 0.1 because of their weak 
lines (see Section |43] |. but this does not affect the result that 
WLQs are a larger fraction of the observed quasar popula- 
tion at higher redshift. The typical Lya+N V EW for nor- 
mal quasars at z > 4.2 is also lower than for lower redshift 
quasars, with median values 46.4 A and 63.8 A, respectively. 
The greater incidence of lower EW objects at higher redshift 
could be caused either by selection effects or by some dif- 
ferential physical process (see below). To investigate the for- 
mer hypothesis, we calculate the selection probability of each 
WLQ based on its redshift, i-band luminosity, spec tral slope, 
and L yg e mission-Une EW usin g the method of iFan et al.] 
(120011) and'Richards et al. I (1200 6a'). Using a distribution func- 
tion of H I absorption, we simulate 200 spectra, compute 
ugriz magnitudes including photometric errors, and apply the 
quasar-selection criteria. The results are shown in the inset 
of Figure |5] We find that the selection probabilities of indi- 
vidual WLQs generally increase with redshift, with average 
probabilities of p = 0.53 at z < 3.5 and /:> = 1.0 at z > 4.5. 
This redshift-dependent selection can be explained by the fact 
that absorption by the Lya forest and Lyman Limit Systems 
(LLSs) equalizes the colors of higher-redshift quasars, in- 
dependent of their emission line properties (see Section [T]!. 
There is also significant dependence on the luminosity and 
spectral slope, as well as emission-line EW, at each redshift 
(e.g.. Fan et al. 2001 , see their Figure 5), but the general result 
is that WLQs are easier to select at higher redshift. The fact 
that the observed WLQ fraction is ^ 6% in samples that are 
unbiased with respect to emission-line strength (e.g., FIRST- 
detected quasars, quasars at z > 4.2) implies that the intrinsic 
WLQ fraction may be significantly higher than the ^ 1 % in- 
ferred from the full SDSS sample. 

In addition to affecting the colors of WLQs, and thus their 
likelihood of selection, absorption by the Lya forest can also 
affect the Lya emission line itself, pushing objects toward 
lower EWs at higher redshift. To have a measurable effect, 
the absorbing material must be located in close proximity to 
the quasar,'^ as is clearly the case for 11 WLQs that have 
PDLAs. However, there is no evidence to suggest that PD- 
LAs, which are generally ass ociated with galaxies in the envi- 
ronrnent of the qua sar je^^g., lM0lleretalJl 19981: lElhson et alj 
I2002t iRussell et al.ll2006h . are more common at higher red- 
shift, and an analysis of the evolution of the intergalactic 
medium (IGM) opacity at z > 3 is beyond the scope of this 
paper. The relevant effect of the Lya forest effect on the 
Lya emission line, thoug h, can be seen in the composite spec- 
trum of z 6 quasars tFan et al.ll2004l) . and it is likely to bias 
the EW measurements of the highest-redshift sources. 

3. IR PHOTOMETRY 

This is the first of three sections (Sections [3] H] |5]l that 
present multiwavelength observations of four z > 4 WLQs: 
SDSSJ1302, SDSSJ1408, SDSSJ1442, and SDSSJ1532. 
These sources were selected early on in the survey to be rep- 
resentative of the handful of WLQs known at the time: their 
Lya emission-line strengths range from undetectable (e.g., 
SDSSJ1532) to weak, but clearly detected (e.g., SDSSJ1442), 

We measure Lya flux at Aiest > 1160 A, so absorbers must be within 
Az = 0.046, corresponding to a peculiar velocity Av = 3100 km s"' or a 
physical distance Z) ~ 8 Mpc at z = 3.5, and Av = 2500 km s"' or D ~ 5 Mpc 
at z = 4.5. 



and two are FIRST radio sources (SDSSJ1408, R = 14; 
SDSSJ1442, R = 26), while the other two have no radio de- 
tections. The FIRST-detected sources also have strong Chan- 
dra X-ray de tections (SDSSJ1408 a„, = -1.54: SDSSJ1442, 
a,„ = -1.42 iShemmer etall I2006h . while SDSSJ1302 and 
SDSSJ153 2 are X-ray weak (a„, = -2.08 and a,„ < -1 .90, re- 
spectivelv: IShemmer et al. ' 2006, 2009). This section presents 
Spit zer mid-IR photometry in the 3.6-24 pm range (Sec- 
tion l3.1b and near-IR photometry in /-, H-, and /iTs-band (Sec- 
tion |32]i, which we use to construct rest-frame 0.1-5 pm 
SEDs of WLQs that can be compared to those of normal 
quasars (Section |33] l. 

3.1. Spitzer Mid-IR Data 

We obtained mid-IR photometry for all four sources as part 
of a Spitzer Space Telescope (Werner et al. 2004) Cycle I Gen- 
eral Observer Program (PID 3221). Observations with the 
IRAC instrument ( Faz io et a l. 2004) were executed in all four 
channels (3.6 /im, 4.5 pm, 5.8 /im, and 8.0 /im) with an inte- 
gration time of 10 00 s per channel. Observations with MIPS 
dRieke et al.ll2004l) were obtained in the 24 /im channel with 
integration times of 1400-1500 s. 

We begin our analysis of the IRAC data using the Basic 
Calibrated Data (BCD) products from version S 14.0.0 of the 
Spitzer Science Center software pipeline. The MOPEX soft- 
ware package was used to construct mosaics from the indi- 
vidual frames, which involves background matching between 
overlapping images and outlier pixel rejection. The native 
IRAC pixel size of 1.2" yields somewhat undersampled data, 
so we take advantage of a subpixel dither pattern to interpolate 
the individual BCD images onto a common grid with a 0.6" 
pixel scale using the Drizzle algorithm. We perform aperture 
photometry on the mosaics using the APER task in IDLPHOT. 
We use 3.6" source apertures with the sky sampled in a 14.4- 
20" annulus. Aperture corrections are calculated based on the 
IRAC point-spread function (PSF) and applied to the flux den- 
sity measurements. The statistical uncertainties in our mea- 
surements are calculated based on Poisson noise in the source 
counts, uncertainty in the background level, and noise from 
sky variations in the source aperture. Except for SDSS J 1442 
in channel 3 and channel 4, the statistical uncertainties for 
all mo saics are smaller than the 5% calibration accuracy of 
IRAC (iReach et al.l2005h . so we quote a 5% uncertainty on all 
of our flux density values. We also calculate array-location- 
dependent photometric corrections for each mosaic; the dif- 
ferences between photometry performed on the corrected and 
uncorrected mosaics are quite small in channels 1 and 2 (1%- 
2 %) and somewhat larger in channels 3 and 4 (l%-5 %). 
We quote the average source flux density from the corrected 
and uncorrected mosaics in Table [3] which is appropriate for 
sources that are flat in through the IRAC bands. We do not 
apply any other color corrections or pixel phase corrections 
because we found them to be < 1 % effects. 

We are unable to perform accurate aperture photometry on 
the IRAC data for SDSSJ1532 because of the presence of a 
star 3.7" away (3.4" east, 1.4" south). Based on the avail- 
able near-IR (Section 13.21 ) and Spitzer data, this source ap- 
pears to be an M star that peaks at the band (/^ ^ 80 /iJy) 
and is fainter at longer wavelengths (/^ ^ 70 piy at IRAC 
channel 1, ^ 45 piy at channel 2). We use the APEX soft- 
ware package to perform point-source fitting photometry on 
the SDSSJ1532 mosaic images (accurate to 5%), and present 
these results in Table [3] 

Our analysis of the MIPS data is based upon the BCD prod- 
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(i2006h) scaled to all of the data points (blue dotted and red dashed lines, respectively). 



ucts from version S 14.4.0 of the software pipeline. We discard 
the first three data collection events (DCENUM=0,1,2) for 
each Astronomical Observation Request to alleviate the first- 
frame effect. We perform additional flat-fielding with flat- 
field.pl in MOPEX to remove a background gradient across 
the frames. We then create mosaics, matching backgrounds 
and removing outlier pixels. We perform aperture photome- 
try on each mosaic using a 6" source radius and a 20-32" sky 
annulus. Aperture corrections are applied to the measured 
flux densities based on the MIPS 24 jivn. PSF. We also per- 
form point-source fitting photometry on all the mosaics, and 
find results consistent with the aperture photometry, within 
the uncertainties. We apply a 4% color correction (i.e., we 
divided the flux density measurements by 0.96), which is ap- 
propriate for a ~ po wer law typical for quasars in this 
rest-fra me frequency range jElvis e t al. 1994i [Richards et alj 
l20061j) . Our results are presented in Table[3] 

3.2. Ground-based Near-IR Data 

We obtained /-band photometry for all four targets and H- 
band photometry for SDSSJ1302 and SDSS J 1408 on 2005 
April 26-27 with the 256 x 256 near-lR camera on the Stew- 
ard Observatory Bok 2.3-m telescope. The camera uses a 
Near-Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) array ( iRieke et al.il 19931) and has a plate scale of 0.6" 
pix"'. The data were taken in nine-position dithered se- 
quences of 60-s exposures. The seeing ranged from 1.3" to 
2.0", and we use 1.8-3.6" source apertures and a 6-12" sky 
annulus for photometry. 



We acquired /Q-band photometry for SDSSJ1302 on 2005 
April 30 and for SDSSJ1408 and SDSSJ1442 on 2005 May 
1 with PANIC ( Martini et al.ll2004l) on the 6.5-m Baade tele- 
scope at Las Campanas Observatory. PANIC utilizes a Rock- 
well 1024 X 1024 IR Hawaii detector with a plate scale of 
0.125" pix"'. The data were taken in five-position dithered 
sequences of 20-s exposures. The seeing ranged from 0.3" 
to 0.5", and we use 0.6-1.9" source apertures and a 2.5-3.1" 
sky annulus for photometry. 

All four objects were observed in H and band, and all 
but SDSSJ1532 were observed in / band on 2005 June 30 
usin g the Near-I R Camera/Fabry-Perot Spectrometer (NIC- 
FPS iHearty et al. 2004) on the Apache Point 3.5-m telescope. 
The camera has a Rockwell Hawaii-IRG 1024 x 1024 CCD 
with a plate scale of 0.273" pix"'. The data were taken in 
five-position dithered sequences in a 20" box surrounding 
the center of the chip, with 60-s exposures in the J band, 
20-s exposures in the H band, and 10-s exposures in the Kg 
band. The seeing during the observations of SDSSJ1302 and 
SDSSJ1408 ranged from 0.7" to 0.9", and the seeing during 
the observations of SDSSJ1442 and SDSSJ1532 ranged from 
1.3 to 1.5". We use 1.1-2.7" source apertures with a 2.7-5.5" 
sky annulus for data obtained in good seeing and a 5.5-6.8" 
annulus for data obtained in poorer seeing. 

For all of the near-lR data, we stack and sky-subtract 
dithered sequences of exposures using procedures included 
in the PANIC IRAF package. For each sequence, we form 
a median sky frame that is scaled to each image and sub- 
tracted. We perform aperture photometry on the science 



8 



Diamond-Stanic et al. 



Rest Wavelength {fim 



Rest Wavelength (/xm) 



10" 



E 

o 

p 10" 

lU 



r SDSSJ1302 




power low (a = 


-0.20) 


blockbody (T= 


870 K) 


combined 










, - ' * 







1 10 
Observed Wavelength {fj,m) 



10" 



10-'- 



r SDSSJ1408 




!x power low (or 
■ •■■^.^ blockbody {T- 


= 0.21) : 


1 150 K) 


^■^.^ combined 

























1 10 
Observed Wavelength (^^m) 



Rest Wavelength (/xm) 



10" 



10" 



SDSSJ1442 



power low (a= 0.02) 
blockbody (T=1240 K) 
combined 



1 10 
Observed Wavelength {fj,m) 



10"' = 



10" 



Rest Wavelength {fin 




: SDSSJ1532 




power low (a = 

blockbody {T= 

combined 


-0.26) 
890 K) 


•• * . 













1 10 
Observed Wovelength (/xm) 



Fig. 7. — Data are the same as Figuie|6] The dotted blue line shows the best-fit power-law model to the short-wavelength data and the dashed red line shows 
the best-fit single-temperature blackbody to the longer wavelength data. The dash-dotted green curve is the sum of the two models. Also included are the spectral 
slope of the power law and the temperature of the blackbody for each source. 



targets and calibrate the measurements using sources in the 
field from the Two Micron All Sky Survey (2MASS) cata- 
log ([Skrutskie et al. 2006), and standard stars from the cata- 
log of Persson et ah (1199 8) that were observed when condi- 
tions were photometric. The uncertainties on these flux den- 
sity measurements are 5%-20%. We quote the photometry 
from each epoch and combine all results into a weighted av- 
erage flux density for each source in each band in Table H) 
Our science targets are too faint to be detected by 2MASS, so 
these are the only photometric epochs available. 

3.3. Rest-frame 0.1-5 i^im SEDs 

The observed SEDs from SDSS / band through MIPS 
24 /im are shown in Figures |6] and |7] Figure |6] also shows 
the SP SS quasar composite spectrum from Vanden Berk et al.l 
(120011) and the mean quas ar SEDs from IeMs et all d 19941) 
and [Richards et al.] (I2006bl) . The data points are from differ- 
ent epochs, so variabiUty may be an issu e, bu t any fluctua- 
tions are likely to be < 10% (see Section l43l l. We normal- 
ize the SDSS composite to match the /- and z-band points, 
and we normalize the template SEDs to provide an overall 
match to all of the data points. The normalized templates are 
reasonably good descriptions of the data for the two radio- 
undetected objects, SDSSJ1302 and SDSSJ1532, in the sense 
that no data point deviates from the template by more than 
30% over the full range of wavelengths from Arest = 1300 A 
to Arest = 5 jim. The two radio-detected objects, SDSS J 1408 
and SDSSJ1442, are brighter at shorter wavelengths by 90% 
and 60% at the / band, and fainter at longer wavelengths 
by 30% and 40% at the MIPS band, respectively, than the 



mean quasar SEDs, but given the ^ 0.5 de x Icr scatter of 
individual quasars aro und the mean SEDs ( Elvis etal.llI9"94t 
[Richards et al.ll20b6al) . these are not significant deviations. 

In Figure |2l we fit the data points with a simple model 
that includes a power-law fit to the short-wavelength (Arest < 
1 /im) flux from the accretion disk (or a relativistic jet), com- 
bined with a single-temperature blackbody fit to the longer 
wavelength thermal emission from hot dust. We allow the 
power-law slope and the blackbody temperature to vary, and 
we fit both simultaneously. We find power-law slopes in 
the range -0.26 < a < 0.21 and dust temperatures in the 
range 870 K < T < 1240 K. This model fits the SEDs of the 
radio-detected quasars better (x^/dof = 1.7 for SDSSJ1408 
and 1.3 for SDSS J 1442) than it does the SEDs of the radio- 
undetected quasars (xV^of = 4.7 for SDSSJ1302 and 5.6 for 
SDSSJ1532). The radio-detected sources are also both bluer 
and have higher temperature dust, which is consistent with 
their redder [5.8] - [8.0] colors. 

4. OPTICAL VARIABILITY 

In this section, we present multiple epochs of op- 
tical spectroscopy (Section 14.1b and /-band photometry 
(Section 14.2b . designed to assess whether SDSSJ1302, 
SDSSJ1408, SDSSJ1442, and SDSSJ1532 exhibit vaiiabihty 
in line or continuum emission that could explain the nature of 
their weak lines. 

4. 1 . Multi-Epoch Spectroscopy 

We obtained spectra for all four sources on 2004 July 7- 
8 with the IMACS spectrograph (.Bigelow & Dressler..2003i) 
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Fig. 8. — Four epochs of optical spectroscopy for SDSSJ1302. The SDSS 
spectrum (R ~ 2000) is binned by a factor of 3, the IMACS spectrum 
(R ~ 3000) is binned by a factor of 5, and the MARS spectrum (R ~ 750) is 
unbinned. Assuming the break associated with the onset of the Lya forest is 
at Arest = 1216 A, we measure z = 4.48 as the systemic redshift. In the second 
panel, we identify an LLS atz = 4.37 with associated Lya absorption, as well 
as weak Si IV and C IV emission lines that are blueshifted by ~ 5000 km s"' 
with respect to ^ = 4.48. In the third panel, we mark the Lya component of a 
strong absorption-line system at z = 4.19, which has metal lines C II A1334, 
Si IV AA1393, 1402, Si II A1526, and C IV AA1548, 1550. There are also two 
weaker absorption-line systems at z = 4.27 and z = 3.98 Lya, Si IV, and C IV. 



on the 6.5-m Baade Telescope at Las Campanas Observatory. 
Observations were made using a 1" slit and a 600 line mm"' 
grism blazed at 7700 A. The spectra were dispersed onto two 
2kx4k chips, one that covered 5690-7660 A and another that 
covered 7720-10090 A. The spectral resolution ranged from 
R - 2500 near A = 6000 A to - 3500 near A = 9000 A. Total 
integration times ranged from 60 minutes for SDSSJ1408 to 
120 minutes for SDSSJ1532. 

Another epoch of spectroscopy was acquired on 2005 May 
8-9 with the MARS spectrograph on the 4-m Mayall tele- 
scope at Kitt Peak. Observations were carried out using a 
1" slit with the 8050-450 grism (450 hnes mm"', blazed at 
8050 A), which afforded wavelength coverage from 5400A 
to 1.1 /im. An OG-550 filter is built into this grating and 
blocks light shortward of A = 5500A. The resolution of the 
spectra is R ^ 750. The total integration time was 40 min- 
utes for SDSSJ1408 and 60 minutes each for SDSSJ1302, 
SDSSJ1442, and SDSSJ1532. 

SDSSJ1302 was also observed on 20 02 Ja nuary 12 with the 
Echelle Spectrograph and Imager (ESI, Sheinis et al .120021) on 
the Keck II telescope for 15 minutes with R ^ 4000. 

We construct flat-field images using dome flats and gener- 
ate wavelength solutions for each object spectrum using lines 
identified from HeNeAr lamp observations. We combine sep- 



FlG. 9. — Three epochs of optical spectroscopy for SDSSJ1408. The reso- 
lution and binning are as in Figure[8] We measure the systemic redshift to be 
z = 3.95, which is lower than the redshift given by SDSS z = 4.01. We iden- 
tify an LLS at z = 3.91 with associated Lya absorption. Weak Si IV and 
C IV emission lines are detected in the second panel; the latter line may 
be blueshifted by ~ 3000 km s"' , but the gap between the IMACS chips 
make this highly uncertain. In the bottom panel, we mark a Lya absorber at 
z = 3.77. 



arate exposures with pixel shifts where appropriate and re- 
move cosmic rays. We generate sensitivity functions using 
standard star spectra and apply the flux calibration to the sci- 
ence targets. One-dimensional spectra were extracted using 
1-2" aperture sizes for the IMACS spectra and from 3^" 
aperture sizes for the MARS spectra, using the best-fit linear 
trace to the two-dimensional spectrum. We make corrections 
for telluric absorption based on extracted standard star spec- 
tra. The calibrated spectra are shown in Figures [SUTTl which 
also include the SDSS spectra for SDSSJ1302, SDSS J 1408, 
and SDSSJ1442 obtained in 2000-2 001, and th e spectra of 
SDSSJ1532 previously published by iFan et alj ([1999) from 
the Double Imaging Spectrograph on the Apache Point 3.5- 
m telescope and the Low-Resolution Imaging Spectrograph 
(LRJS) on the Keck II telescope. We record the spectral slope 
and the EW of Lya, as well as Si IV and C IV in cases where 
the lines are detected, for all spectra in Table |5] 

4.2. Multi-Epoch Photometry 

SDSS /-band photometry (iFukugita et alJ 119961) was ob- 
tained for a ll four sources on 2005 May 11 with OPTIC 
( iTonrv et al . 2004) on the 3. 5-mWIYN telesc ope and on 2006 
April 21-22 with 90Prime (Wifliams et al. 2004) on the Stew- 
ard Observatory 2.3-m Bok telescope. OPTIC consists of two 
2k X 2k CCDs wifli aplate scale of 0.14" pix"', while 90Prime 
has four 4kx4k CCDs with a plate scale of 0.45" pix"'. The 
seeing during the observations ranged from 1.2" to 2.0". 
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Fig. 10. — Three epochs of optical spectroscopy for SDSSJ1442. The res- 
olution and binning are as in Figure [8] We measure the systemic redshift to 
be z = 4.51, and we identify an LLS at the same redshift. The rest-frame 
LyQ EW increases with time from 9 A in the SDSS spectrum to 16 A in the 
MARS spectrum. We also detect Si IV emission in this latter epoch. In the 
bottom panel, we mark a Lya absorber at z = 4.47. 



Fig. 11. — Four epochs of optical spectroscopy for SDSSJ1532. The spec- 
tra in the top two panels are taken directly from Fan et al. 1 1999), and the 
resolution and binning of the IMACS and MARS spectra are as in Figure [8] 
A redshift of z = 4.62 is consistent with both the onset of the Lya forest and 
the location of the LLS. In the second panel, we mark az = 4.58 system with 
Lya and Ly/3 absorption. The rest-fraine EW of Lya is less than 3 A in all 
epochs. 



For the OPTIC data, we construct a night sky flat using a 
combination of the science images, rejecting the five highest 
and five lowest values for each pixel, and averaging the re- 
maining values. For the 90Prime data, we construct a flat-field 
image using 10 twilight sky flats, rejecting the two highest and 
one lowest values of each pixel and averaging the remaining 
seven pixel values. The flat-field images are scaled to their 
median value before we apply them to the data. 

We perform aperture photometry on each science target and 
on nearby SDSS stars in the field that are (1) brighter than 
the science target and (2) faint enough to be in the regime 
where the CCD response is linear to better than 1 % (< 30, 000 
ADU). These "standard" stars are all fainter than / = 17, 
have SDSS photometry that is accurate to ^ 1% ( Ivezi c et"an 
l200l . and are within 5' of the science target. We use 2.8^.2" 
source apertures with a 4.9-5.6" sky annulus for the OPTIC 
data and 2.7-5.4" source apertures with a 6.8-9.0" sky annu- 
lus for the 90Prime data. We measure the flux of the science 
target, with appropriate aperture corrections, relative to each 
standard star, and then calculate the weighted average of these 
measurements, which we quote in Table|6]and use to construct 
light curves in Figure [12] We also include SDSS photometry 
in our analysis, which was obtained in 1999-2001 (see Ta- 
ble©. 

4.3. Variability, Weak Broad Lines, and Redshifts 

Over the seven-year baseline of these observations, we find 
only modest photometric and spectral variations. As can be 
seen in Figure [121 the /-band continuum emission fluctua- 



tions are less than ±10% for all sources. Similarly, there is 
not much variability in the Lya, Si IV, and C IV emission 
lines, which remain weak in every spectrum in Figures Isl-ITTI 
The most notable variation can be seen in Figure [TO] where 
SDSS J 1442 has a clear detection of Lya in all epochs, as well 
as a detection of Si IV in the MARS spectrum. The rest-frame 
EW of Lya increases with time, ranging from 9 A in the SDSS 
spectrum to 16 A in the MARS spectrum, placing it outside 
our WLQ definition (EW < 10 A). There are also detections 
of weak Si IV and C IV emission lines in the Keck spectrum 
of SDSSJ1302 (see Figure |8]l and in the IMACS spectrum of 
SDSSJ1408 (see Figure|9]l. These lines are not detected in the 
other available spectra, which have lower S/N. 

The lack of strong emission lines makes the task of measur- 
ing accurate redshifts more difficult for WLQs than for nor- 
mal quasars. The primary spectral feature in these sources 
is the break at the onset of the Lya forest, and we estimate 
the redshift for each WLQ by assuming that this feature is 
at Ai-est = 1216 A. These estimates are consistent with pre- 
viously published redshifts for SDSSJ1442 (z = 4.51) and 
SDSSJ1532 (z = 4.62), both of which also have an LLS at the 
same redshift (see Figures [10] and fTTT i. but they are different 
from the redshifts in the Fifth Data Release Quasar Catalog 
for SDSSJ1302 (we find z = 4.48 compared with z = 4.4684) 
and SDSSJ1408 (we find z = 3.95 compared with z = 4.0075). 
The LLSs in the latter two spectra indicate somewhat lower 
redshifts (z = 4.37 and z = 3.91, respectively), but this is not 
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Fig. 12. — Light curves from i'-band photometry spanning a baseline of at least six years for each WLQ. The dotted line in each panel is the weighted average 
of all the data points, which corresponds to i = 19.84 for SDSSJ1302, i = 18.91 for SDSSJ1408, i = 20.03 for SDSSJ1442, and / = 19.80 for SDSSJ1532. The 
synthesized /-band magnitudes from the spectroscopy epochs are also consistent with the mean values, but they have larger uncertainties than the photometry and 
we do not plot them for clarity. The fluctuations around the mean flux density are less than ±10% for all sources. 



necessarily a problem because one does not always expect an 
LLS to be exactly at the redshift of the quasar 

However, the weak emission lines detected in SDSSJ1302 
(second panel of Figure[8]) are consistent with a lower redshift 
if we interpret them as Si IV and C IV. The central wave- 
lengths of these lines correspond to z = 4.38 and z = 4.39 or 
blueshifted velocities v = 5500 km s"' and v = 4900 km s"' rel- 
ative to z = 4.48. We also tentatively measure a C IV blueshift 
in SDSSJ1408, but the IMACS chip gap makes the cen- 
tral wavelength of C IV highly uncertain. The C IV emis- 
sion line is typically b lueshifted wi th respect to the systemic 
quasar redshift (e.g.. iGaskelll 11982 ). but the amplitude of 
the velocity is < 1000 km s~' in most cases. For example, 
[Richards et al.l (l2002bl) measured the velocity shift of C IV rel- 
ative to Mg 11 for ~ 800 quas ars and found 824 ±511 km s"' 
(mean and dispersion). Shen et all (120081) found similar re- 
sults, including very few C IV blu eshifts > 3000 km s~' , 
in a sample of 15,000 quasars. [Richards et al.l (l2002bl) 
also found that the sources with larger C IV blueshifts have 
smaller EWs, so the interpretation that SDSSJ1302 has high- 
ionization emission lines blueshifted by ^ 5000 km s"' is ex- 
treme, but perhaps feasible given the weakness of its lines. 
Alternatively, the systemic redshift could be lower than our 
estimate of z = 4.48, implying that the line and continuum 
emission have been strongly absorbed at Arest = 1216 A, which 
is the case for some PDLAs with Zabs — Zem (see Section|2|. 



5. POLARIZATION AND RADIO CONTINUUM 



In this section we present optical polarimetry (Sec- 
tionEB and VLA observations (SectionESli of SDSSJ1302, 
SDSSJ1408, SDSSJ1442, and SDSSJ1532 designed to assess 
whether they could be powered by synchrotron emission from 
a relativistic jet, as is the case for BL Lacs. 

5.1. Multi-Epoch Optical Polarimetry 

We obtained optical polarization measurements on 2005 
May 13-16 and 2006 May 4-5 using the CCD Imag- 
ing/Spectropolarimeter (SPOL. [Schmidt et alj|1992ah on the 
Steward Observatory 2.3-m Bok Telescope. The instrument 
was used in imaging polarimetry mode, which provides a 
51" X 51" field of view with ^ 0.5" pixels. In this config- 
uration, the key optical elements are a rotatable A/2 wave- 
plate and a Wollaston prism, which passes both orthogonal 
polarizations and splits them into separate beams that are fo- 
cused onto a 1200 X 800 SITe CCD detector. We obtained 
data using sequences of exposures that produce two images 
for both beams, which can be analyzed to measure the Q and 
U Stokes parameters. We calibrated the data by (1) inserting 
a NICOL prism into the beam, (2) observing an unpolarized 
standard star, a nd (3) measuring polarized standard stars from 
Schmidt et at] ( 11992b ). 

In 2005 May, we used an / filter; in 2006 May, we used 
an R6% filter, which is designed to transmit light redward 
of 6800 A and has > 90% transmittance beyond 7200 A. 
The seeing during the observations ranged from 1.5" to 2.5". 
Additional imaging polarimetry data for SDSSJ1532 were 
obtained on 1999 May 18 (these data were published by 
[Fan et al. .l999,) and 2000 January 8. These observations used 
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an R12 filter designed to transmit light redward of 7200 A. 

We trim and bias-subtract the images using standard rou- 
tines in IRAF. We create a master flat-field image for both the 
Q and U sequences by median combining the output of three 
Q dome-flat sequences and three U dome-flat sequences. Af- 
ter flat fielding, we perform arithmetic operations on images 
from both beams to produce Q and U images, as well an / 
Stokes image corresponding to the total flux from the source. 
We perform aperture photometry on these images to measure 
the normalized Stokes parameters q = Q/I and u = U /I. We 
used 3" apertures and a 5-10" sky annulus for all the images. 
Absolute flux calibrations are not necessary for our polariza- 
tion measurements, so we make no effort to correct for source 
flux extending outside the photometric aperture. 

We calculate the linear polarization Pobs = \/ (f- + u?- and the 
polarization angle d = 0.5 arctan(M/g') for each science tar- 
get. The results are shown in Table |7] and discussed in Sec- 
tion 15.31 We quote both the observed polarization Pobs and 
the polarization with a first-order correction for statistical bias 
dWardle & K ronberg^ l974 . P = (P^s " f^?) the latter value 
is used in our analysis. 

5.2. VLA Radio Continuum Data 

We obtained observations with the VLA in the A-array con- 
figuration (Thomps on et aLl il980) on 2006 May 10. All four 
sources were observed in the L band (A ~ 20 cm) with the LI 
frequency combination, which consists of 2 x 50 MHz inter- 
mediate frequencies (IFs) centered ?Av= 1365 MHz and v = 
1435 MHz. To measure radio spectral slopes, the two sources 
with FIRST radio detections (SDSSJ1408, SDSSJ1442) were 
also observed in C-band (A ^ 6 cm) with 2 x 50 MHz IFs cen- 
tered at 1/ = 4835 MHz and v = 4885 MHz . There were 22 
antennas in operation during the observations, and the longest 
antenna baseline was ^ 36 km, which provided an angular 
resolution of 1 .4" FWHM in the L band and 0.4" in the C 
band. The total on-source integration times were 150 minutes 
for SDSSJ1532, 100 minutes for SDSSJ1302, and 25 min- 
utes at both L and C bands for SDSSJ1408 and SDSSJ1442, 
respectively. Primary and secondary calibrator sources were 
observed periodically, and absolute calibration was tied to ob- 
servations of 3C 286. The data for our science targets were 
calibrated using standard routines in AlPS. Several antenna 
baselines were affected by RFl noise in the L band, and we 
flagged the UV data with anomalous amplitudes. The data 
were cleaned and maps were made with the IMAGR task 
within AlPS. The morphologies of the two detected sources 
(SDSSJ1408, SDSSJ1442) are consistent with point sources. 
The flux from each was measured using the IMEAN task with 
circular apertures at least twice the radius of the beam. The 
resultin g flu x densities are given in Tableland discussed in 
SectionlSJl 

5.3. Optical Polarization, Radio Loudness, and Spectral 

Slope 

We do not detect any significant optical polarization from 
SDSSJ1302 in either 2005 May or 2006 May. The la uncer- 
tainty in this measurement is 0.68%, and we set a 3cr upper 
limit of P < 2.0% for this WLQ. We do detect a low level of 
polarization at both epochs for SDSSJ1408 and SDSSJ1442, 
with weighted average values of P = 1.50% ±0.37% and 
P = 2.43% ± 0.63%, respectively. From measurements of 
foreground stars, we estimate the value of interstellar polar- 
ization along the fine of sight to SDSSJ1408 and SDSSJ1442 



to be no greater than ^0.1%. We detect polarization from 
SDSSJ1532 in all four epochs of observation, although never 
with > 3(7 significance in a single epoch. This source is at 
a fairly low Galactic latitude, and we measure interstellar 
polarization P 1%. Taking the weighted average of the 
Q and U values from all four epochs of observations for 
SDSSJ1532, we find P = 1.72% ±0.44%. If we group 
the two older observations (1999 May and 2000 January) 
separately from the two newer observations (2005 May and 
2006 May), we find 2.91% ±0.91% and 1.30% ±0.50% 
respectively. Three of the WLQs in our full z > 3 sample 
(SDSSJl 14153. 34H-021924.3, SDSSJ121221.56H-534127.9, 
and SDSS J 12374 3.08-H630144.8) were observed by 
ISmith et alj (l2007h . who found P = 0.90% ± 0.64%, 
P = 1 .32% ± 0.58%, and P = 1 .21 % ± 0.78% respectively. All 
of these values are below the nominal level f or the definition 
of highly polarized objects (P > 3%, e.g., Impey & Tapial 
,1990) and within the range for optically s elected quasars 
witho ut a synchrotron component (e.g., iBerriman et alj 
IT990I) . 

We do not detect SDSSJ1302 or SDSSJ1532 in our L-band 
radio observations. The maps for these two sources have Icr 
noise levels of 35 pJy beam"' and 28 piy beam"', respec- 
tively. We calculate limits on the radio-loudness parameters 
for these two sources using flux densities at rest-frame 6 cm 
determined from the 3fT upper limits at observed-frame 20 cm, 
assuming a,- = -0.8, and flux densities at rest-frame 2500 A 
extrapolated between the observed frame 1 .2 and 1 .6 /im (J 
and H bands) flux densities. We find R < 2.5 for SDSSJ1302 
and < 1.8 for SDSSJ1532. A flatter assumed radio spec- 
tral slope would place even tighter limits of R (e.g., for ar = 
the limits drop to P < 1.7 and P < 1.2, respectively). We 
measure flux densities for SDSSJ1408 at the L band and C 
band of 0.89 ± 0.07 mJy and 0.46 ± 0.06 mJy, correspond- 
ing to a radio spectral slope ar = -0.52. The L-band mea- 
surement is 25% or 1.9(7 smaller than the FlRST-integrated 
flux density (1.18 mJy, rms noise 0.14 mJy beam"'), so there 
may be moderate variability, but it is not highly significant. 
The flux densities for SDSSJ1442 at the L band and C band 
are 0.83 ±0.07 mJy and 0.45 ±0.06 mJy, corresponding to 
a,- = -0.49. The L-band measurement is 56% or 6.3cr smaller 
than the FIRST integrated flux density (1.87 mJy, rms noise 
0. 15 mJy beam"'), which indicates significant variability. The 
values of the radio spectral slope are both near the division be- 
tween steep-spectrum and flat-spectrum radio sources at ar = 
-0.5. We calculate the radio-loudness parameters for these 
two sources using their 20 cm flux densities, their measured 
radio spectral slope, and an extrapolation between the 1 .2 and 
1.6 ^im flux densities. We find R = 10.0 for SDSSJ1408 and 
P = 29.3 for SDSSJ1442. 

6. THE NATURE OF WEAK EMISSION-LINE QUASARS 

As shown in Section|2l the line luminosities of high-redshift 
WLQs are 4x fainter than those of normal z > 3 quasars, 
and their continuum luminosities are 40% brighter There are 
two hypotheses that are consistent with this result: (1) WLQs 
are intrinsically less luminous than normal quasars in terms 
of both line and continuum emission, but a relativistic jet 
beamed toward us amplifies their continua; or (2) WLQs have 
the same intrinsic continuum properties as normal quasars, 
but some physical process, either a lack of line-emitting gas 
or obscuration along the line of sight, causes the observed 
Lya and other UV emission lines to be weak. In this section. 
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we discuss how our results fit within these two hypotheses, 
and what physical processes may be at work. 

6.1. Arguments Against Continuum Boosting 

6.1.1. UV-IR Properties 

We find that all four WLQs with Spitzer photometry show 
emission from hot (T ^ 1000 K) dust. This places strong 
constraints on any continuum boosting because relativistic 
jets have no effect on thermal dust emission. We can thus 
rule out continuum boosting for the radio-undetected WLQs 
(SDSSJ1302, SDSSJ1532), whose UV-IR SEDs closely 
match those of typical quasars in all respects, including their 
ratio of power law to thermal dust emission. The radio- 
detected sources (SDSSJ1408, SDSSJ1442) have weaker dust 
emission by factors of 1.5-1.9 (in the MIPS band) relative to 
the mean quasar SEDs, but this is well within the factor of 
2-3 scatter of normal type 1 quasar SEDs. Thus their UV- 
optical emission may be boosted by a factor of ^ 2, but no 
more, which is not sufficient to explain the extreme weakness 
of their lines 

Further evidence against continuum boosting comes from 
the fact that WLQs do not exhibit strong optical polariza- 
tion. The low levels of polarization observed in SDSSJ1408, 
SDSSJ1442, and SDSSJ1532 are probably intrinsic, but are 
too small to imply that this polarization comes from syn- 
chrotron emission, as is the case for BL Lacs. Radio-selected 
BL Lacs are found to be highly polarized (P > 3%) ^ 90% 
of the time (i.e., in a single epoch, one would expect that 
9/10 radio-selected BL Lacs would show P > 3%) and X-ray- 
selected BL Lacs are found to be polarized roughly half of the 
time (e.g., Jannuzi et al. 1994). The fact that no WLQ is found 
to be highly polarized in several epochs of observations indi- 
cates that these objects are significantly less polarized than 
even X-ray-selected BL L acs and consistent wi th the polar- 
izations of normal quasars dBerriman et alJ[T990l) . 

There is also no evidence of strong optical variability in 
WLQs."' The fluctuations that are seen in Figure [T2]are_con- 
sistent with those of normal quasars (e.g., Vanden Ber k et alJ 
12004). We do see variability in both Lya EW and radio flux in 
SDSSJ1442, but the lack of any corresponding optical contin- 
uum fluctuations in this source argues against the continuum 
boosting scenario. There is no evidence to suggest that the 
increase in its line strength with time is a reverberation effect, 
although our temporal sampling is not ideal since the light 
crossing time of the broad-line region (BLR) is expect ed to 
be ^ 5 yr in the observe d frame at this luminosity (e.g.. lLaorl 
ll998l:lBentzetalJl2009h . The drop in the radio flux is signifi- 
cant, and it is not required that the radio and optical emission 
of a jet vary in a synchronized manner, but BL Lacs are found 
to be variable at all wavelengths on a variety of timescales, so 
the lack of any significant optical continuum variability indi- 
cates that this emission is not likely coming from a jet. 

6.1.2. Radio Properties 

If the physical process causing weak line emission were 
continuum boosting by relativistic jets, we would expect a 

IStalin & SriananJ <2005l) measured a ~ 20% drop in flux in R and / 
bands for SDSSJ1532 between 2000 and 2001 and claimed more dramatic 
fading relative to the SDSS i'-band photometry in 1999. Their results are 
called into question by the second epoch of SDSS photometry in 2001 (see 
Table|6), which shows only a 10% drop relative to 1999. Their claim of a 
~ 50% flux difference between the 1999 SDSS i'-band epoch and their 2000 
R-band photometry is erroneous because they fail to account for the fact that 
the i?-band flux is heavily absorbed by the Lya forest. 



large fraction of WLQs to be radio loud. Instead, we find 
no statistical difference between the distribution of radio- 
loudness parameters for WLQs and normal quasars. The 
precise number of sources with R > 10, a value which is 
often used to describe radio loudness'^, is not well con- 
strained for either normal quasars or WLQs because most 
have radio upper limits that are still consistent with R > 10. 
The SDSS and FIRST data do definitively indicate, how- 
ever, that 92% of WLQs and somewhere between 93.3% and 
96.1% of normal quasars have R < 40, in contrast to BL 
Lacs, the majority of which have R > 100 (see, e.g.. Fig- 
ure 5 of IShemmer et al.1 l2009l) . Almost ah WLQs are ei- 
ther radio-quiet or radio-moderate, indicating that if their con- 
tinua are boosted, the effect is monochromatic (i.e., roughly 
equal in the radio and rest-frame UV) and distinct from 
the je t mechanism at work in BL Lacs. Shemm er et alJ 
(2009) discussed the potential as sociation of WLQs with 
radio - weak BL Lacs (e.g., Londish et al.)l2004l;ICollin2e et alJ 
120051; lAnderson et"ai. 2007; Plotkin et al. 2008) and pointed 
out that the lack of typical (i.e., radio-loud) BL Lacs at 
high redshift makes it difficult to connect the two phenom- 
ena. Further evidence that we are not seeing pole-on ra- 
dio jets in WLQs comes from the ^ -0.5 radio spectral 
slopes for SDSS J 1408 and SDSS J 1442, which are signifi- 
cantl y steeper than the t ypical slopes for BL Lacs, a,. ~ 0.3 
(e.g.. 1Stickel etan[T99Th . 

6.2. The Remaining Possibilities 

A variation on the continuum boosting hypothesis in- 
volves gravitational lensing, where WLQs could either be (1) 
strongly lensed galaxies or (2) normal quasars whose con- 
tinuum emission has been microlense d by a star in an inter- 
vening galaxy. She mmer et alJ (l2006h ruled out the strongly 
lensed galaxy hypothesis for WLQs with strong X-ray detec- 
tions on the basis of their X-ray-to-optical flux ratios, which 
are typical for quasars. We additionally rule it out for all 
four WLQs in Figure |6] because their UV-IR SEDs match 
those of normal quasars. The microlensing hypothesis has 
received some attention in the literature, and several authors 
have used the variability and EW distributions of quasars 
to put constraints on the properties of lensing objects (e.g ., 
iDalcanton et"ai][T994l; IZackrisson et al.ll2003l;IWiegertll2003l) . 
The characteristic timescale for mi crolensing is ^ 10 yr for a 
stellar lens in a foreground galaxy ( lGouldlfl995h . so we can- 
not rule out microlensing for WLQs, but there is no evidence 
of fading continua over 6-7 yr of observations. 

The variety of arguments against continuum boosting as 
a cause for the weak emission-line strength of WLQs, cou- 
pled with evidence against lensing, implies that WLQs are 
a rare, unique population at high redshift. However, there 
are several objects at lower redshift whose physical prop- 
erties may be related. One of these is the z = 0.94 radio- 
quiet quasar PG 1407-H265, whose properties are described 
in detail by iMcDowell et al.l (119950 . This object has very 
weak Lyg emission (EW= 8 A ) and does not exhibit po- 
larization (Berriman et al. 1990") or optical-UV variability. 
Blundell et al. (2003) and Gallo (2006) observed a factor 
of ~ 2 variability in the radio and X-ray bands, simil ar to 
the radio variability we see in SDSSJ1442 (Section l53b , 
but this seems to be unrelated to its optical-UV continuum 
flux. It also has a weak Mg 11 emission line (EW = 24 A), 

" A more stringent definition is R > 100, in which case objects with 10 < 
R < 100 are refen'ed to as radio-moderate. 
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a somewhat stronger Ha line (EW =126 A), and unusu- 
ally strong Fe II emission. Interestingly, its weakly detected 
C IV line (EW = 4 A) is blueshifted by - 4000 km s"' 
with respect to Lya, similar to the blueshift for SDSSJ1302 
discussed in Section 14.31 Other objects we are aware of 
with weak Lya emission and blues hifted C IV line s are 
SDSS.T152156.48+52023 8.4 (z = 2.2 I.Tust et al.1 l2007h and 
HE 0141-3932 (z = 1.8, Reimers et alJI2005h . IRichards et al.1 
(|2002b) argued that observed C IV blueshifts are due to a lack 
of flux in the red wing of the emission line, and their discus- 
sion of this phenomenon in the context of cloud-based and 
accretion-disk-wind models for the BLR is pertinent here. 

Under the hypothesis that WLQs have normal quasar con- 
tinuum properties, but that some physical process causes the 
emission lines to be weak, several of the possible interpreta- 
tions for PG1407H-265 mentioned by McDowell et al. ( 1995]) 
are relevant for WLQs: (1) the BLR could have anomalous 
properties or a low covering factor, (2) an exceptional geom- 
etry could cause the BLR to see a continuum that is different 
from the one that we see, or (3) the BLR could be covered by 
a patchy BAL region that does not affect the continuum. In 
the context of explaining physical processes that would result 
in weak or absent broad lines, N icastro et al.l (1200 3) discussed 
a scenario where the BLR forms via accretion disk instabili- 
ties at a critical radius where gas pressure begins to domi- 
nate over radiation pressure; at low accretion rates / luminosi- 
ties the BLR would move toward smaller radii and eventually 
cease to exist. Similarly, Laor (2003) pointed out that if the 
BLR cannot survive at line widths Av > 25,000 km s"', a 
minimum luminosity / accretion rate is implied, and he dis- 
cussed how the outer and inner boundaries of the BLR may 
be set by suppression of line emission by dust and thermal 
processes. Czernv et al. (2004) went a step further and cal- 
culated the minimum radius, minimum Eddington ratio, and 
maximum line width for which the BLR is expected to exist 
in ADAF and disk-evaporation models. However, the associ- 
ation of high-redshift WLQs with low accretion rates, which 
would explain their weak lines in the above theoretical sce- 
narios, is ruled out by the fact that their continuum properties 
(X-ray, UV, optical, IR, radio) are comparable to those of nor- 
mal quasars. 

Another relevant object at low redshift (z = 0.192) is the 
radio-quiet quasar PHL 1811, a narrow-line Se yfert 1 galaxy 
(NLSl) with weak Lya and C IV emission dLeighlv et al.l 
|2007b), weak X -ray e mission, and a steep X -ray spectrum 
(aav = -2.3, r = 2.3, iLeighlv et al.1 12007al) . iLeighlvet al.1 
(l2007b ) showed that the weak high-ionization emission lines 
in this quasar can be explained by its soft SED in the sense that 
a lack of high-energy photons prevents typical gas photoion- 
ization processes from occurring. They speculated that such 
weak X-ray and UV line emission may be associated with a 
high accretion rate, which is often invoked as a physical in- 
terpretation for NLSls (e.g.. Boil er et al .1119961) . Similarly, 
the local NLSl NGC4051 (z = 0.002) exhibited weak X-ray 
emission and a weak He II A4686 broad emission line (ion- 
ization potential 54.4 eV) during the final months of a three- 
year monitoring campaign by Peterson et al. (2000), while 
its broad H/3 line remained strong. The authors speculated 
that the inner part of the accretion disk may have become 
advection-dominated during this period, suppressing the X- 
ray and far-UV continuum and also the high-ionization lines. 

We do not have constraints on the strength of the low- 
ionization lines (e.g., Mg II, Balmer lines) in WLQs, but 



further observations are warranted to test where these lines 
are stronger than the high-ionization lines, as is the case for 
PG 1407H-265, PHL 1811, and NGC 4051 (in its low X-ray 
flux state). In models that describ e the BLR in terms of a 
disk wind (e.g.. lMurrav et al.] |1995). the low-ionization lines 
are produced in the accretion disk and the high-ionization 
lines are produced in the outflowing wind, so the UV emis- 
sion lines could be suppressed either by an abnormal pho- 
toionizing continuum or by a process that prevents the disk 
wind itself from forming; either of these could be associated 
with a high accretion rate. Evidence of the accretion rate be- 
ing the physica l driver of t he Baldwin e f fect is presented by 
iBaskin & Lao3 ( [2004) and lBachev etal.1 ( 12004) . but it is not 
clear if such an inverse relationship between t he C iv EW and 
accretion rate persists at higher luminosities. IShemmer et al.l 
(1^09) explored the possibility that WLQs could be extreme 
quasars with high accretion rates, and correspondingly steep 
X-ray spectra, by jointly fitting the available X-ray data; they 
found a slope that is consistent with those of normal radio- 
quiet quasars, but higher-quality X-ray spectra are required to 
test this hypothesis properly. 

Finally, it is worth considering the effects of absorption on 
the observed LyaH-N V EWs of WLQs. As discussed in Sec- 
tion |2] there is evidence of strong intervening absorption in a 
fraction of the WLQ sample. While we flag those WLQs with 
obvious PDLA systems, absorption by lower column densi- 
ties of material could also affect the remainder of the sample. 
In Figure [3] the decrease in the Lya/N V ratio and the red- 
ward migration of the peak of the Lycu-N V feature as one 
moves toward lower EWs could be explained by absorption 
that affects not just the blue side of Lya, but also the peak 
of the line and emission redward of the peak. Such behav- 
ior is also seen, although to a lesser ex tent, in high-redshif t 
quasar composite spectra presented by Dietrich et al.l (120021) 
and Fan et al. (2004). There is no evidence of correspond- 
ing C IV absorption, however, so the nature of the absorption 
would have to be different than in B ALs and systems wit h 
intrinsic narrow-line absorbers (e.g.. ICrenshaw et"an 120031) . 
One scenario that would explain the absorption of Lya and not 
metal li nes would be the infall of pristine gas from the IGM 
(e.g., Barkana & Loeb 2003). There is certainly evidence of 
strong IGM H I opacity blueward of Lya, and we cannot rule 
out absorption at Arest > 1216 A in some WLQs, but it does 
not explain the C IV weakness, and we conclude that most 
WLQs are likely to have intrinsically weak emission lines. 

7. CONCLUSIONS 

We have identified a sample of 74 high-redshift (z > 3) 
SDSS quasars with weak UV emission lines and obtained IR, 
optical, and radio observations for four such WLQs at z > 4. 
We find that the optical continuum properties of WLQs are 
similar to those of normal quasars, while their Lya + N V line 
luminosities are significantly weaker, by a factor of ~ 4 on av- 
erage. The radio-loudness parameters of WLQs are also sta- 
tistically indistinguishable from those of normal quasars. The 
strong spectral break caused by Lya forest absorption makes 
it easier to optically select WLQs at higher redshift, and the 
high WLQ fraction at z > 4.2 6% of quasars) indicates that 
they may be quite common in general, although difficult to se- 
lect in most quasar surveys. All WLQs with IR observations 
show evidence of hot (T 1000 K) dust emission, which rules 
out hypotheses that invoke relativistic continuum boosting to 
explain their weak line emission. Their Arest = 0. 1-5 fj,m SEDs 
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closely resemble the mean quasar SED. They are not highly 
variable or polarized in the optical, and their weak emission 
lines persist through several epochs of spectroscopy. In ad- 
dition, the radio spectral slopes of radio-detected WLQs are 
also significantly steeper than those of BL Lacs. 

This evidence supports the interpretation that WLQs have 
the same intrinsic continuum properties as normal quasars, 
but that some physical process associated with the line- 
emitting gas results in weak lines. We are not able to dis- 
criminate among various possibilities that could cause the ex- 
ceptional emission-line weakness in WLQs, but they include 
a low BLR covering factor, an abnormal geometry that ex- 
ploits the difference between our line of sight to the accretion 
disk and the BLR, or suppression of high-ionization emission 
lines perhaps related to a high accretion rate. An important 
observational test of the latter scenario will require near-lR 
spectroscopy to search for the low-ionization emission line 
Mg II, as well as Bal mer lines, Fe II, and [O Ill1 to facili- 
tate comparison to the iBoroson & GreenI (119921) Eigenvector 
1 concept. Further study is justified as WLQs offer insight 
into emission-line physics and our overall understanding of 
the quasar phenomenon. 
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TABLE 1 

Parameters for all z > 3 Quasars 



NAME (SDSS J) 


Lya+N V" 


f7 


C IV 


fj 








uniform** 


BAL' 


000046.41+011420.8 


87.2 


3.2 


71.4 


2.6 


-2.08 


3.76 


< 18 





2 


000051.56 + 001202.5 


23.8 


1.3 






-1.14 


3.88 


66 





2 


000135.51-004206.7 


58.3 


1.2 






0.74 


3.58 


< 17 





2 


000221.11+002149.3 


23.9 


0.5 


16.5 


1.5 


-0.21 


3.06 


132 








000238.41-101149.8 


117.3 


2.9 






-0.96 


3.94 


< 28 








000252.72-000331.0 


62.7 


4.7 






-1.02 


3.68 


67 








000300.35 + 160027.6 


63.7 


2.1 


15.2 


2.4 


-1.18 


3.68 




1 





000303.35-105150.6 


89.0 


1.2 


79.8 


2.1 


0.03 


3.65 


< 11 





3 


000316.38-000732.4 


60.8 


1.9 


30.9 


1.7 


-0.08 


3.18 


< 29 








000335.20+144743.6 


27.6 


2.5 






-0.17 


3.48 




1 


1 



Note. — TablelTlis published in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its form and content. 
^ EW [A]. 

^ Calculated as described in Section [2. 1.2| The upper limits use the FIRST catalog detection limit, a 5<t threshold that includes CLEAN bias. 

Uniform selection flag in Fifth Data Release Quasar Catalog: (0) not selected as a primary quasar target by the final target selection algorithm given bv lRichards erarH20Q2^ , (1) 
selected as a primary quasar target by the final target selection algorithm. 

^ BAL flag. (0) not identified as a BAL in either lTVump et all 120061) or[Gibson et al. 1 2009) catalog. (1) identified as a BAL in lTrump et alj f2006|) catalog only, (2) identified as a 
BAL in lGibson et aU 1 .2009.) catalog only, (3) identified as a BAL in both lTrump et aU 1 .2006.) and . Gibson et aU 1 .2009.) catalogs. 
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TABLE 2 

Sample of Weak Emission-Line Quasars 



JNAMh (SUSS J) 


EW (A)"* 


z 


ob 
K 


unifo 


004054.65-091526.7 


4.5 


4.98 


< 16 





00542 L42-01092L6 


12.3 


5.09 


< 24 





010802.90-010946.1 


12.3 


3.37 


< 14 





022337.76 + 003230.6 


8.6 


3.07 


57 





025646.56 + 003858.3 


12.8 


3.47 


< 12 





080523.32 + 214921. l''''= 


5.6 


3.46 


< 10 




082059.35 + 561022.0'''' 


12.2 


3.64 


<23 




083122.57+404623.3 


8.8 


4.88 


< 11 




084434.15 + 224305.2''-'= 


11.7 


3.12 


< 13 




091738.90 + 082053.9'' 


11.2 


3.25 


< 13 




093306.88 + 332556.6 


11.9 


4.56 


< 16 




095 108.76 + 3 14705.8 


5.8 


3.03 


127 




1 n9QAQ 80 4- ^in^V? 1 8'^ 
luzy^y.ou + DUJ / J L.o 


14 


'^ 1 Q 


^ 79 
^ ZZ 




inn/in ^^^^nioin od 




T C7 


^19 

< Iz 




1 0360 1 12 + 084948 4 


9.4 


3.68 


< 15 




104650.29 + 295206.8 


6.3 


4.27 


< 14 




105049.28 +441144.8^ 


5.3 


4.32 


< 13 




113354.89+022420.9'* 


15.0 


3.99 


<7 




113415.21 + 392826.1 


13.8 


4.83 


<29 




11 3729.42 +375224.2« 


14.7 


4.17 


59 




114153.34+021924.3 


3.3 


3.48 


17 




114412.77 + 315800.8 


7.0 


3.23 


< 11 




114434.60 + 510317.8 


9.8 


3.91 


< 24 




114958.54 + 375115.0 


13.3 


4.31 


<34 




115254.97 + 150707.7'' 


12.8 


3.33 


<6 




115308.45 + 374232.1'' 


15.2 


3.03 


<29 




115906.52+133737.7'' 


11.9 


3.98 


< 2 




115933.53+054141.6'' 


10.5 


3.29 


< 20 




120059.68 + 4009 13.1"= 


6.9 


3.37 


64 




121221.56 + 534127.9 


8.2 


3.10 


< 5 




121812.39 + 444544.5 


12.6 


4.52 


< 19 




1 2202 1 .39 + 092 1 35 .8''''= 


12.9 


4.11 


< 11 




1 9?^SQ IS -t- 1 1 7800 1 ti'^ 


1 1 s 


4 12 


^ 1 8 




19744i 77 + '575921 3 


15.1 


4.30 


< 18 




1231 16.08 + 41 1337.3 


10.7 


3.84 


< 18 




123132.37 + 013814.0 


6.3 


3.23 


84 




123315.94 + 313218.4 


12.9 


3.22 


106 




1 9'^szl^ 1 Q -1- 1 9'^^l7^ 7^' 


10 5 


^ 9 1 


<' 7S 

^ ZJ 






4 6 


^ 47 

J. HZ 


^ 7 




1 ZM-ZvJ^.Z / + DZ J / iZ. i 


1 9 9 
IZ.Z 


^ ^7 


78 
zo 




iZJJvJo. / J + i JUoU4.V 


13.9 


3.63 


< 7 




lzDjiy.lU + 4j4iJz.cS 


11.3 


3.53 


< 19 






/.O 


A An 
4,4/ 


< ZJ 




1 J) W J) J Z . f z U Z 1 7 UU . J 


in Q 


4. DO 


^ ZJ 








1 Q 1 
J. 61 


< is 






10.1 


3 77 


<; 7 




133146.20+483826.3^ 


15.3 


3.74 


< 11 




133422.63 + 475033.5 


11.6 


4.95 


< 16 




134453.51+294519.6 


9.1 


4.71 


< 14 




134521.39 + 281822.2 


6.3 


4.08 


< 24 




140300.22 + 432805.3 


14.5 


4.70 


< 21 




140850.91+020522.7 


7.0 


4.01 


14 




141209.96 + 062406.8 


2.4 


4.47 


771 




1413 1 8.86 +450522.9''''' 


3.4 


3.11 


1782 




142103.83 + 343332.0'= 


7.0 


4.91 


< 11 




1 47 1 44 QS + 3i n 1 i 4''.'= 


3.1 


4.56 


< 9 




1A99S7 f^/i -1- '^7'^S^7 

iHZZ J / .UO T J / JoU / .D 


s 9 
o.z 


3 16 


^ 7 1 

^ Zi 






7.0 


3 72 


< 17 






19^ 


J. 1 V 






L't'tZ J 1 . / 1 iJ 1 i 1.7 J J .J 


1 1.4 


4 5 1 


25 




150220.46 + 465233.5 


11.4 


4.26 


< 18 




150739.67-010911.0 


12.9 


3.14 


< 15 




152200.14+413741.7 


9.3 


3.24 


17 




155203.30 + 352440.4 


4.6 


3.04 


<8 




155645.31 + 380752.7 


14.9 


3.32 


< 13 




160336.64+350824.2 


5.3 


4.46 


<6 




161122.44 +414409.6'''^ 


8.6 


3.13 


<7 




163411.82+215325.0 


6.7 


4.53 


< 10 




210216.52+104906.6'' 


14.7 


4.18 






214753.29 - 073031.3 


14.5 


3.15 


18 




223827.17 + 135432.6'''= 


1.2 


3.52 






225246.43 + 142525.8 


11.1 


4.90 
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TABLE 2 — Continued 



NAME (SDSS J) 


EW (A)" 




uniform'^ 


233255.71 + 141916.4 


10.9 


4.75 


1 


233446.40-090812.2'' 


12.6 


3.32 


102 



EW(Lya+N V) 

^ Calculated as described in Section [2.1.2| The upper limits use the FIRST catalog detection limit, a 5cr threshold that includes CLEAN bias. 

Uniform selection flag in Fifth Data Release Quasar Catalog: (0) not selected as a primary quasar target by the final target selection algorithm given bv lRichards etaj] I2002;j) . (1) 
selected as a primary quasar target by the final target selection algorithm. 

EW(C IV) > 10 A. 
^ In lProchaska et al] I20Q8I) catalog of Proximate Damped Lya Absorbers. 



TABLE 3 
Spitzer PHOTOMETRY 



Name 


3.6 [^Jy] 


fj 


4.5 [MJy] 


cr 


5.8 [^l}y] 


cr 


8.0 [^l}y] 


cr 


24.0 [^Jy] 


cr 


SDSSJ1302 


78.3 


3.9 


62.6 


3.1 


62.9 


3.1 


84.6 


4.2 


470 


24 


SDSSJ1408 


86.4 


4.3 


76.3 


3.8 


80.5 


4.0 


140.6 


7.0 


575 


29 


SDSSJ1442 


37.9 


1.9 


33.0 


1.7 


35.2 


2.7 


67.8 


3.7 


191 


24 


SDSSJ1532 


83.5 


4.2 


70.8 


3.5 


72.1 


3.6 


99.8 


5.0 


534 


27 



TABLE 4 
Near-IR Photometry 



Name 


Telescope/Instnjment 


Date 


J [/i Jy] 


cr 




cr 


K [/iJy] 


cr 


SDSSJ1302 


Bok2.3m/256x256 


2005 Apr 26 


59.4 


12.1 


67.7 


13.1 






SDSSJ1302 


Magellan/PANIC 


2005 Apr 30 










59.7 


6.0 


SDSSJ1302 


APO 3.5m/NICFPS 


2005 Jun 30 


54.3 


7.0 


69.6 


10.2 


58.6 


8.6 


SDSSJ1302 


Weighted average 




56.3 


9.3 


68.9 


11.4 


59.1 


7.8 


SDSSJ1408 


Bok2.3m/256x256 


2005 Apr 27 


105.3 


9.9 


101.5 


17.1 






SDSSJ1408 


Magellan/PANIC 


2005 May 1 










92.0 


5.2 


SDSSJ1408 


APO 3.5m/NICFPS 


2005 Jun 30 


113.4 


8.5 


94.3 


10.6 


92.0 


6.9 


SDSSJ1408 


Weighted average 




111.3 


9.4 


95.1 


11.6 


92.0 


6.8 


SDSSJ1442 


Bok2.3m/256x256 


2005 Apr 27 


38.2 


7.4 










SDSSJ1442 


Magellan/PANIC 


2005 May 1 










34.4 


2.9 


SDSSJ1442 


APO 3.5m/NICFPS 


2005 Jun 30 


32.1 


5.3 


38.6 


6.4 


34.7 


5.1 


SDSSJ1442 


Weighted average 




34.2 


6.3 


38.6 


6.4 


34.7 


4.1 


SDSSJ1532 


Bok2.3m/256x256 


2005 Apr 27 


64.0 


8.3 










SDSSJ1532 


APO 3.5m/NICFPS 


2005 Jun 30 






80.6 


6.7 


66.1 


7.3 


SDSSJ1532 


Weighted average 




64.0 


8.3 


80.6 


6.7 


66.1 


7.3 



High-Redshift Weak Emission-Line Quasars 
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TABLE 5 
Spectroscopic Parameters 



Object 


Instrument 


Date 


a 


Lya 




Si IV 


C IV 








EWrest [A] EWrest [A] 




SDSSJ1302 


SDSS 


2001 Mar 26 


-0.45 


2.5 








SDSSJ1302 


Keck/ESI 


2002 Jan 12 


-0.61 


2.4 




2.5 


4.7 


SDSSJ1302 


IMACS 


2004 Jul 8 


_A AQ 


Z.J 








SDSSJ1302 


MARS 


2005 May 9 


-0.49 


2.3 








SDSSJ1408 


SDSS 


2001 Mar 25 


-0.57 


7.0 












OfM\d. Till 7 
ZUU^ Jul / 


-0.55 


2.0 




7.4 


o.u 


SDSSJ1408 


MARS 


2005 May 9 


-0.55 


6.1 








SDSSJ1442 


SDSS 


2000 Apr 29 


-0.23 


9.1 








SDSSJ1442 


IMACS 


2004 Jul 7 


-0.33 


12.9 








SDSSJ1442 


MARS 


2005 May 9 


-0.15 


15.7 




7.4 




SDSSJ1532 


APO 


1999 Mar 13 


-1.50 


1.1 








aUbaJ 1532 


KbCK/LKlS 


1999 May 13 


-0.96 


1.2 








SDSSJ1532 


IMACS 


2004 Jul 8 


-0.72 


2.5 








cr\ccTi CIO 


MARS 


2005 May 8 


-1.00 


1.1 












TABLE 6 














SDSS i-BAND Photometry 












Instrument 


Date 




a [/uJy] 






SDSSJ1302 


SDSS 


1999 Mar 21 




41.8 


1.1 






SDSS Jl 302 


WIYN/OPTIC 


2005 May 11 




42.4 


1.3 






SDS5iJ13()2 


Bok/90Prime 


2006 Apr 21 




41.2 


2.1 






SDSSJ1408 


SDSS 


2000 May 4 




104.5 


1.8 






SDSSJ1408 


WIYN/OPTIC 


2005 May 1 1 




91.7 


2.2 






SDSSJ1408 


Bok/90Prime 


2006 Apr 21 




98.3 


4.0 






SDSSJ1442 


SDSS 


1999 Mar 20 




38.3 


1.2 






SDSSJ1442 


SDSS 


2000 May 4 




35.1 


1.1 






SDSSJ1442 


WIYN/OPTIC 


2005 May 11 




35.3 


1.9 






SDSSJ1442 


Bok/90Prime 


2006 Apr 21 




32.5 


1.1 






SDSSJ1532 


SDSS 


1999 Mar 21 




44.1 


1.0 






SDSSJ1532 


SDSS 


2001 Jun 15 




39.4 


1.2 






SDSSJ1532 


WIYN/OPTIC 


2005 May 11 




43.7 


1.4 






SDSSJ1532 


Bok/90Prime 


2006 Apr 21 




46.7 


1.1 
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Diamond-Stanic et al. 



TABLE? 
Polarization 



Name 


Date 


q 




u 






pb 


crp 






SDSSJ1302 


2005 May 13 


0.67 


1.09 


-0.08 


1.07 


0.68 




1.09 


176.7 




SDSSJ1302 


2006 May 4 


0.11 


0.96 


-0.49 


0.77 


0.50 




0.78 


141.0 




SDSSJ1302 


0.36 


0.72 


-0.35 


0.63 


0.50 


- 


0.68 


157.8 




SDSSJ1408 


2005 May 13 


2.05 


0.84 


-0.19 


0.72 


2.06 


1.88 


0.84 


177.4 


n.i 


SDSSJ1408 


2006 May 4 


0.87 


0.41 


0.83 


0.38 


1.21 


1.14 


0.40 


21.9 


10.0 


SDSSJ1408 ^ 


0.00 


0.12 


-0.04 


0.11 


0.04 




0.11 


132.6 




SDSSJ1408 f 




1.42 


0.38 


0.59 


0.34 


1.54 


1.50 


0.37 


89.7 


5.7 


SDSSJ1442 


2005 May 13 


1.02 


1.26 


1.85 


1.20 


2.11 


1.73 


1.21 


30.6 


20.1 


SDSSJ1442 


2006 May 4 


1.33 


0.75 


2.19 


0.72 


2.57 


2.46 


0.73 


29.3 


8.5 


SDSSJ1442 <= 


-0.08 


0.01 


0.11 


0.01 


0.13 


0.13 


0.01 


62.5 


2.8 


SDSSJ1442 f 




1.54 


0.64 


1.99 


0.62 


2.51 


2.43 


0.63 


26.1 


7.4 


SDSSJ1532 


1999 May 18 


1.45 


1.13 


2.23 


1.05 


2.66 


2.44 


1.08 


28.5 


12.7 


SDSSJ1532 


2000 Jan 8 


-1.33 


1.65 


3.48 


1.75 


3.73 


3.30 


1.74 


55.4 


15.1 


805511532 = 




-0.83 


0.01 


0.00 


0.01 


0.83 


0.83 


0.01 


89.9 


0.4 


SD55J1532 ' 




1.37 


0.93 


2.72 


0.90 


3.05 


2.91 


0.91 


31.6 


9.0 


5D55J1532 


2005 May 13 


0.46 


1.35 


3.88 


1.55 


3.90 


3.59 


1.54 


41.6 


12.4 


SD55J1532 


2006 May 4 


-0.09 


0.49 


1.08 


0.53 


1.08 


0.94 


0.53 


47.4 


16.1 


SDSSJ1532' 


-0.99 


0.13 


0.39 


0.14 


1.06 


1.06 


0.13 


79.3 


3.7 


SDSSJ1532 f 




0.95 


0.48 


1.01 


0.52 


1.39 


1.30 


0.50 


23.4 


11.1 


SDSSJ1532 8 




1.04 


0.43 


1.44 


0.45 


1.78 


1.72 


0.44 


27.1 


7.4 



Note. — The data in columns 3-9 are all percentages. 

6 = 0.5 arctan(M/^) 

(79 =28.65 (ffp/P) 

^ Measurement of interstellar polarization, 
f 

Weighted average corrected for interstellar polarization. 
^ Weighted average of four epochs corrected for interstellar polarization. 



TABLE 8 
VLA Continuum Observations 



Name 


Lband 
[mJy] 


rms noise 
[mjy/beam] 


Cband 
[mJy] 


rms noise 
[mJy/beam] 


a 


R 


SDSSJ1302 




0.035 








< 2.5 


SDSSJ1408 


0.891 


0.074 


0.461 


0.060 


-0.52 


10.1 


SDSSJ1442 


0.830 


0.066 


0.446 


0.060 


-0.49 


29.3 


SDSSJ1532 




0.028 








< 1.8 



